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The electromagnetic structure ol the lightest hadrons, proton, pion, and kaon, is studied by 
high precision measurements of their form factors for the highest timelike momentum transfers of 
IQ^I = s = 14.2 and 17.4 GeV^ Data taken with the CLEO-c detector at ^ = 3.772 GeV and 
4.170 GeV, with integrated luminosities of 805 pb~^ and 586 pb~^, respectively, have been used to 
study e^e~ annihilations into 7r"^7r~, K'^K~ , and pp. The dimensional counting rule prediction that 
at large the quantity Q^F{Q^) for pseudoscalar mesons is nearly constant, and should vary only 
weakly as the strong coupling constant, as(Q^), is confirmed for both pious and kaons. However, the 
measurements are in strong quantitative disagreement with the predictions of the existing QCD- 
based models. For protons, it is found that the timelike form factors continue to remain nearly 
twice as large as the corresponding spacelike form factors measured in electron elastic scattering, in 
significant violation of the expectation of their equality at large Q^. Further, in contrast to pious 
and kaons, a significant difference is observed between the values of the corresponding quantity 
\Q'^\GM{\Q^\)/fJ-p for protons at |Q^| = 14.2 GeV^ and 17.4 GeV^. The results suggest the constancy 
of IQ^IGiv/dQ^DZ/ip, instead, at these large \Q^\. 



PACS numbers: 13.40Gp,14.20Dh,14.40Be,14.40Df 

Knowledge of the quark-gluon structure of the only 
stable baryon, the proton, and the lightest mesons, the 
pion and kaon, is of great interest for both nuclear and 
particle physics. Important questions about the size of 
the proton, the composition of its spin, and the large 
difference between its spacelike and timelike form fac- 
tors remain open. Timelike form factors of pions and 
kaons, which are needed for the precision determination 
of the hadronic loop contribution to the muon g — 2 
anomaly [l!, are poorly known. Spacelike form fac- 
tors of pions and kaons needed for the understanding of 
nuclear and hypernuclear forces are difficult to measure 
at large momentum transfers, and can only be obtained 
by analytic continuation of timelike form factors 3] . To 
meet these needs, precision measurements of timelike 
form factors at the highest possible momentum transfers 
are needed. In this Letter we report measurements of 
the form factors of pions, kaons, and protons with much 
higher precision, and for much larger timelikc momentum 
transfers than before [1, Q . 

Earlier measurements of proton form factors for large 
timelike momentum transfers {Q^ < 0) made by the Fer- 
milab E760/E835 pp — )■ e~^e~ experiments for jQ^j = 
8.84 - 13.11 GeV^ 0], and the CLEO e+e" pp 
measurements at |Q^| = 13.48 GeV^ 0], revealed that 
the timelike form factors are nearly twice as large as 
the corresponding spacelike form factors, a result in 
strong disagreement with the expectation of their equal- 
ity [6| at asymptotically large \Q'^\- The measurement 
of pion and kaon timelike form factors by CLEO at 
= 13.48 GeV^ Q revealed that while the dimensional 
counting ruleprediction of a as/|Q^| variation of the 
form factors [6| was apparently confirmed, the measured 



form factors were factors 4 — 8 larger than predicted. Fur- 
ther, the ratio Fy^dQ'^l)/ Fk{\Q'^\) was also found to be 
nearly twice as large as the QCD prediction that it should 
be equal to the ratio of the squares of the decay constants, 
f^/ fx at large These large differences from the- 

oretical predictions raise important questions about how 
valid the asymptotic predictions are at the measured mo- 
mentum transfers, and make it imperative to extend the 
measurements to larger momentum transfers. 

We use data taken with the CLEO-c detector, which 
has been described in detail before (sj, to study the re- 
actions e^e~ — >■ pp, 7r"'"7r~, and K~^K~ . The main body 
of the data comprise of e^e^ annihilations at ccnter-of- 
mass energies Js = 3772 MeV {\Q^\ ^ 14.2 GeV^) and 
4170 MeV (|Q2| = 17.4 GeV^), with integrated lumi- 
nosities of £ = 805 pb~^ and L — 586 pb"^, respec- 
tively. Data from a mini-scan at the average ^fa = 
4010.4 MeV ((IQ^j) = I6.O8 GeV^, £ = 20.7 pb'^) and 
yfs = 4260 MeV (IQ^j = 18.25 GeV^, £ = 12.9 pb'^) 
have also been analyzed. 

No measurements of the branching fractions for non- 
D~D two-body decays of either ^(3772) or ?/;(4160) to 
light hadrons exist [9|, but they can be estimated us- 
ing the perturbative QCD (pQCD) prediction that the 
hadronic decays of iplnS) states scale with the princi- 
ple quantum number n in the same way as the dilep- 
ton decays. The measured values, S(^(3770) e+e~) 
and ^(-0(4160) — > e+e"), are ~ 10^ times smaller than 
B{ip{2S) e+e") This leads to the estimate that 
the branching fractions for tt+tt", K^K~ , and pp decays 
of -0(3770) and ■0(4160) are approximately 0.9 x 10"*^, 
9 X 10~*, and 3.2 x 10^*", respectively, based on the mea- 
surements for V'(2S') decays 0, [13] ■ This leads to the 
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estimates that with more than 5.2 milhon ■0(3772) and 
■0(4160) formed in the present measurements, the num- 
bers of resonantly produced tt+tt^, K^K^, and pp are 
~ 0.04, 0.4, and 1.8, respectively. In other words, the 
resonance contribution to their yield is expected to be 
vanishingly small, and the counts observed in the present 
analysis can be entirely attributed to the timelike elec- 
tromagnetic form factors. 

The event selection and particle identification for the 
analysis reported in this Letter are similar to those in 
our earlier publication Q. The reconstructed events are 
required to have two charged tracks, zero net charge, and 
I cos 6*1 < 0.80. Each charged particle is required to satisfy 
the standard criteria for track quality and origin of the 
track at the interaction point. In order to develop data- 
independent particle identification criteria, and to deter- 
mine event selection efficiencies, Monte Carlo (MC) sam- 
ples were generated for e+e" — > h'^h^ , h — p,n, K using 



the EvtGen generator [ll|, and e+e 



1+1-, I 



using the Babayaga generator [121. For e 



and K^K~, the Monte Carlo samples were generated 
with sin^ 9 angular distributions, since for pseudoscalar 
mesons, the differential cross sections are 



daois,9)„Jdn = iay8s)l3l\FUsW 



(1) 



where a is the fine structure constant, c/3„i is the meson 
velocity in the laboratory frame, and Fm{s) is the form 
factor for timelike momentum transfer at s = \Q'^\. For 
e+e" — pp the Monte Carlo samples were generated with 
both (1 + cos^ 9) and sin^ 6 angular distributions, since 

daois, ey^/dn = {a^/As)Pp[\GM{s)\^{l + cos^ 9) 

+ t\Ge{s)\^ sinH], (2) 

where Ge{s) and Gm{s) are the electric and magnetic 
form factors of the proton, respectively, and r = Arrip/s. 

Monte Carlo generated momentum distributions for 
different pairs of hadrons (h^h^) and leptons (l^l^) 
show that the peaks for individual particles are broad- 
ened by detector resolution, and they overlap and develop 
large tails at low momenta, mainly due to FSR. Further, 
the several orders of magnitude larger QED yields of lep- 
tons (e^e~ and fi^fi~) overwhelm the hadron peaks. 
The first step of analysis therefore consisted of mini- 
mizing the lepton contributions and the radiative tails. 
This preselection of events was done by removing events 
with E'cc (energy loss in the central calorimeter)/p(track 
momentum)> 0.85 (which reduced e^e~ by a factor 
~ 10^), removing those with any hits in the muon cham- 
bers (which reduced muons by a factor ^ 10^), and re- 
quiring that the vector sum of the momenta of the pair 
be J2p+,- < 60 MeV/c (which greatly reduced the ra- 
diative tails). To develop further event selection criteria, 
we use the variable Xh = {E{h^) + E{h~'))/^, as in 
Ref. i. 



For particle identification in the RICH detector, 
we define the parameter Li, which is the likelihood 
that a particle corresponds with a given hypothesis 
of being of species, i, based on the Cherenkov pho- 
tons detected in the RICH detector. We combine 
this with dE/dx information from the drift chambers, 
X^{dE/dx), defined as x'^ (dE/dx) = [{dE/dx) 

measured 

(d£'/da;)oxpcctod]^/CT^oasurod' to coustruct the joint vari- 
able AC{i,j) to distinguish particle type i from contam- 
inant particle type j. 



A£(i,j) = -2[log(L,)-log(ij)] 



+ [x\dE/dx),-x\dE/dx)j]. (3) 



Monte Carlo simulations show that the default values 
of AC{i,j) < are very effective in distinguishing be- 
tween the desired particle i and the contaminant particle 
j. Accordingly, we require AC{p, /i), AC{p, e), AC{K, /i), 
AL{K,e), AC{K,7r), AC{K,p), A/:(7r,e), AL{t:,K), 
AC{'r:,p) < 0. Because pions and muons have similar 
masses, a AC{i,j) cut is not effective in distinguishing 
between pions and muons. It was shown in Ref. ^i] that 
Ecc can be very effective in distinguishing muons, which 
suffer only ionization loss in the CC crystals, from pions, 
which suffer additional energy loss due to hadronic in- 
teractions in the crystals. We therefore impose the addi- 
tional requirement of Ecc > 350 MeV. With this addi- 
tional requirement, muon contamination in the pion peak 
is found to be less than 1%, and the Xt^ k.p distributions 
for both = 3772 MeV and 4170 MeV are aU free of 
contaminants, as shown in Fig. 1. 

Since the peaks in Fig. 1 have essentially no back- 
grounds, we obtain the number of counts A'^(s) and 
Nk{s) as counts in the range X^^^k — 1.000 ± 0.015, 
and Np{s) as counts in the range Xp ~ 1.000 ± 0.010. 
We note that the Xp distribution for = 3772 MeV 
in Fig. 1(e) has a definite tail in the low Xp region. As 
shown by the unshaded histogram in Fig. 1(e), this arises 
from pp from the decay of i/'(2S', 3686) produced by ini- 
tial state radiation (ISR), and is clearly observed when 
the net vector momentum in the events is increased to 
< 150 MeV/c. 

The observed cross sections are obtained as cro(s) = 
N{s)/[e{s) ■ >C(s)], where e(s) is the MC-determined effi- 
ciency, and C{s) is the integrated e^e^ luminosity. The 
efficiencies at ^/s = 3772 MeV and 4170 MeV are respec- 
tively 16.2% and 16.2% for tt+tt", 60.2% and 60.9% for 
K+K-, and 71.3% and 68.7% for pp. These cross sec- 
tions are corrected for ISR to obtain the Born cross sec- 
tions using the method of Bonneau and Martin [l3] • At 
y/s — 3772 MeV and 4170 MeV the correction factors are 
respectively 0.797 and 0.796 for tt+tt", 0.817 and 0.809 
for K+K-, and 0.806 and 0.800 for pp. The Born cross 
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TABLE I. Cross sections for e+e 



K+K , and pp 



for e+e annihilations at ^ = 3772 MeV and 4170 MeV, and 
the corresponding form factors of pion, kaon, and proton. 



FIG. 1. Distributions of = [E{h+) + E{h-)]/y/s ior h = n, 
K, p for = 3772 MeV and 4170 MeV data. The vertical 
dashed hnes bracket the Xh region in which counts are ac- 
cepted. The open histogram in panel (e) corresponds to pp 
from tp{2S) ISR excitation (see text). Arrows indicate ex- 
pected positions of contaminants. 



sections are related to the angle integrated cross sections 
ob{s)p = (47raV3s)/3p[|GA/(s)P + {t I2)\Ge{s)\^] 



(4) 
(5) 



The cross sections as and the form factors F^,_r- and Gm 
(assuming Ge = Gm) derived from them are listed in 
Table L The systematic uncertainties listed in the table 
are as described later. In the Table, and elsewhere in this 
Letter, the errors indicated in the parentheses are in the 
corresponding last digits of the main values. 
Pion Form Factors — Timelike form factors of pions, 
measured by e+e~ — >■ 7r-'"7r~, have been reported by 
Babar for 0.09 GcV^ < \Q^\ < 8.7 GeV^ by the ISR 
method At the largest |Q^|, these results have up to 
±100% errors, and were based on < 10 counts/bin. Our 
earlier measurement [3 at = 13.48 GeV^ had only 
26(5) observed counts, and resulted in F^r (13.48 GeV^) — 
0.075(9). The present resuhs, F^(14.2 GeV^) = 0.065(2) 



TT+TT 
pp 


,K+K- N^,K 
Np 


as (pb) WF^,K 
as (pb) IO'Gm 


\Q'^\Gm / fJ'p 


Vs = 


3772 MeV, |Q^j = 


14.2 GeV^ 




K+K 
pp 


661(26) 6.36(25)(36) 0.65(1)(2) 
1564(40) 3.95(10)(22) 0.54(1)(1) 
213(15) 0.46(3)(3) 0.88(3)(2) 


0.92(2)(3) 
0.76(1)(2) 
0.64(2) (2) 


V~s = 


4170 MeV, \Q''\ = 


17.4 GeV^ 




11+ TT 

K+K 
pp 


218(12) 
644(25) 
92(10) 


2.89(16)(16) 0.48(1)(1) 
2.23(9)(12) 0.44(1)(1) 
0.29(3)(2) 0.76(4)(2) 


0.84(2) (2) 
0.77(2)(2) 
0.82(4)(2) 



and F^{17A GeV^) = 0.048(1) are based on 661(26) 
and 218(12) observed counts, respectively. These are 
listed in Table I. In Fig. 2 we plot our results together 
with all previous results, including the indirect result for 
IQ^i = M^{J/'ijj) [l^. As shown in the figure, all mea- 
surements with IQ^I > 9 GeV^ are in excellent agree- 
ment with the dimensional counting rule prediction of a 
l/IQ^I variation of the form factors at large momentum 
transfers. Also shown in the figure are three illustrative 
theoretical predictions. The behavior of the predic- 
tion of the QCD sum rule inspired model [16] disagrees 
strongly with the data. The pQCD prediction of Gous- 
set and Pire [l3| is nearly a factor two smaller than our 
measurements, and the latest AdS/QCD prediction by 
Brodsky and de Teramond [l3| reproduces the data be- 
low 5 GeV^, but falls to 2/3 of the observed values for 
IQ^I > 5 GeV^. Czyz et al. [l9| have shown that the mea- 
sured F^^dQ^I) at IQ^I > 5 GeV^ can be parameterized 
in a VDM approach, but only by including hypothesized 
radial ps, pi, resonances. 

Kaon Form Factors — The results of our present 
measurements at IQ^I = 14.2 GeV^ and 17.4 GeV^ are 
listed in Table I. In Fig. 2, we show our present re- 
sults for kaon form factors along with early results for 
\Q^\ < 4.4 GeV2 ^2^], the indirect result for \Q^\ = 
M'^{J/'ip) [2l|, and our previous measurement at \Q'^\ = 
13.48 GeV2]5|. As for pions, aU results for IQ^I > 9 GeV^ 
follow the predicted as/\Q'^ \ behavior of the form factors. 
No theoretical predictions for kaon timelike form factors 
exist. An empirical fit to the data has however been made 
by Czyz et al. 19], but it requires an "infinite tower" of 
hypothetical p, u, and <j) resonances. 

Proton Form Factors — In Fig. 2, we show the pre- 
viously published proton form factors for timelike mo- 
mentum transfers. At the largest \Q'^\ all these mea- 
surements were severely limited by statistics, with 
and 1 counts at \Q'^\ = 14.4 GeV^ [i], 16(5) counts at 
]Q2| = 13.48 GeV^ and 8(5) counts in the region 
jg2| = 13.0 - 14.4 GeV^ We observe 213(15) and 

92(10) counts at \Q^\ = 14.2 GeV^ and 17.4 GeV^, re- 
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FIG. 2. Summary of form factor results. The solid points in the pion and proton panels are from BaBar ISR measurements [14l . 
[2^ . the open triangles are from FNAL pp measurements [J, the open circles are from the CLEO measurements the open 
squares at |Q^| = 16.1 GeV^ and 18.3 GeV^ are from our small statistics data sets, and the solid squares are from the present 
measurements at — 14.2 GeV'^ and 17.4 GeV^. The theoretical curves for pions are from References [16,17,18]. The solid 
curves illustrate the arbitrarily normalized variation of as for tt and K, and q| variation for protons. 



spectively. This allows us to examine, for the first time, 
the IQ^I dependence of the proton form factors sensi- 
tively. As listed in Table I, and shown in Fig. 2, we find 
that in disagreement with the dimensional counting rule 
prediction of the weak Og variation of \Q'^\G m{\Q'^\) / l^-p, 
its value 0.64(3) GeV^ at \Q^\ = 14.2 GeV^ is 22(4)% 
smaller than 0.82(5) GeV^ at 17.4 GeV^. We note, how- 
ever, that |Q^|Gm(|Q^|)/Mp is essentially constant, be- 
ing = 4.5(2) X 10^2 Qgv2 and 4.7(2) x IQ-^ QeV^ for 
^ 14.2 GeV^ and 17.4 GeV^, respectively. This is 
reminiscent of the fact that for the spacelike momentum 
transfers for the equivalent form factors, _F'i(Dirac) and 
F2(Pauli), it is QF2/F1 which is found to be constant, 
rather than Q'^F^/Fl as predicted 0- 

Recent polarization measurements of the spacelike 
form factors of proton at JLab have revealed that 
Ge{Q'^)/Gm{Q'^) monotonically decreases with increas- 
ing [23j. We fit the summed differential cross sec- 
tions we measure at ^ 3772 MeV and 4170 MeV 
with Eq. 2, and obtain Gm = (0.85 ± 0.07) x 10"^ 
Ge = (0.71 ± 1.17) X 10-2, and Ge/Gm = 0.83+i];^^ 
at the average (IQ^I) = 16.1 GeV^. This shows that we 
have little sensitivity to Ge, and even with our larger 
statistics we can not determine Ge/Gm- The results 
listed in Table 1 and shown in Fig. 2 are assuming 
GeHQ^D = GmUQ^D- If Gb(|Q2|) ^ is assumed, 
Gm(IQ^I) would increase by - 6% at \Q^\ = 14.2 GeV^, 
and by - 5% at \Q^\ = 17 A GeV^. 

Systematic Uncertainties — As in our previous publi- 
cation 0, we estimate uncertainties of 1% in trigger effi- 
ciencies, 2% in tracking efficiencies, 1% in luminosity, and 
0.2% in radiative corrections. In addition, we estimate 
the total uncertainty due to variation of '^p, Ecc, and 
AC to be < 5% for tt, K, and p at both ^ = 3772 MeV 
and 4170 MeV. This brings the total systematic uncer- 
tainty to 6%. 

To summarize, we have made form factor measure- 



ments for pions, kaons, and protons for the highest 
timelike momentum transfers of \Q'^\ = 14.2 GeV^ and 
17.4 GeV^ with nearly five times higher precision than 
prior measurements in the literature. For pions and kaons 
we find that the dimensional counting rule prediction of 
as/Q^ variation of the form factors with jQ^j holds very 
well. However, the existing theoretical predictions for 
pions fail by large factors to predict the magnitude of 
the form factors. We find F^{\Q^\)/ Fk{\Q^\) = 1.21(3) 
and 1.09(4) for \Q^\ = 14.2 GeV^ and 17.4 GeV^, respec- 
tively. These are in agreement with .F^dQ^DZ-F/f (IQ'^I) = 
1.19(17) measured at 13.48 GeV^ [5], and in large dis- 
agreement with the asymptotic prediction that they 
should be equal to the ratio of the squares of the decay 
constants, f^/fj^ = 0.67(1). 

For protons, we find that the timelike form factors con- 
tinue to be a factor two or more larger than the cor- 
responding spacelike form factors, as shown in Fig. 2. 
We find the unexpected result that \Q'^\Gm{\Q'^\)/ f^p at 
|Q2| = 14.2 GeV^ is 22(4)% smaller than at \Q^\ = 
17.4 GeV^. The difference is suggestive of the near- 
constancy of |Q2|Ga/(|Q2|)/^p, instead. 

The overall conclusion of the present investigation is 
that the asymptotic predictions of QCD-based models 
are not realized even at momentum transfer as large as 
18 GeV^, and that theoretical understanding of timelike 
form factors of hadrons is still lacking at the quantita- 
tive level, and our precision measurements provide new 
challenges for theory. 

This investigation was done using CLEO data, and as 
members of the former CLEO Collaboration we thank it 
for this privilege. This research was supported by the 
U.S. Department of Energy. 
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